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Abstract 


Using data collected with the CLEO II detector at the Cornell Electron 
Storage Ring, we present new measurements of the branching fractions for 
D'^ —> KsK^ and iL^Tr"*". These results are combined with other 

CLEO measurements to extract the ratios of isospin amplitudes and phase 
shifts for D KK and D —> iCvr. 
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Strong final-state interactions (FSI) in nonleptonic weak decays of hadrons obscure the under¬ 
lying weak interactions. The problem is particularly acute for the D meson, as its mass lies in 
a resonance-rich region M- Elastic (i.e. nir stays as vrvr) and inelastic FSI rotate the isospin 
amplitudes 0. These isospin amplitudes may be inferred by combining measurements of branch¬ 
ing fractions. This Letter reports new measurements of the D~^ KgK^ and 
branching fractions. We combine these results with previous CLEO measurements of branching 
fractions [||-|^ to obtain the first measurement of the isospin amplitudes and phase shift difference 
for D ^ KK and improved values of these quantities for D —> LCtt. 

The CLEO H detector is designed to measure charged particles and photons with high 
efficiency and precision. This analysis is based on 3.12 fb~^ of data collected at the T(45) resonance 
and 1.72 fb~^ 60 MeV below the T(4S'). Hadronic events are selected by requiring at least three 
charged tracks, a total detected energy of at least 0.15 Ec.m., and a primary vertex within 5 cm 
along the beam (z) axis of the interaction point. 

Candidate Ks mesons are detected in the Ks mode. They are reconstructed by 

combining pairs of oppositely charged tracks, each with an impact parameter in r — iji of greater 
than four times the measurement uncertainty. The track pair must also pass a cut based on the 
the distance in z between the two tracks at their r — cj) intersection point. The invariant mass of 
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the track pair must be within 15 MeV of the known Ks mass. 

Charged pion and kaon candidates must pass minimum track-quality requirements. To reduce 
combinatoric background in the D~^ —> channel, we require that the specific ionization 

{(lE/dx) of the candidate be within 3 standard deviations {a) of that expected for a kaon. 
Tighter cuts are applied on the candidates in the —> KsK~^ mode because of a large 

background from D~^ —> Kst^~^ decays. The measured dE/dx must be within 2a of that expected 
for a kaon and at least 0.25a lower than that expected for a pion. 

We then reconstruct candidates from the Ks, K~^, and vr"’' candidates in the signal modes 
D~^ —> KsK^ and Ks'x'^, and the normalization mode In the Kst^~^ 

mode, we observe a large background from events in which a Ks candidate is combined with a 
random slow pion. We therefore require cos{9ks) < 0-8, where 9ks is the angle between the Ks in 
the rest frame and the D~^ direction in the laboratory frame. This requirement is also imposed 
on the D~^ —> KsK~^ mode. 

We require that every candidate also be a product of the decay D*^ —> . The low- 

momentum provides a clean tag for the D'^. Pairs of electromagnetic showers detected by 
CLEO’s CsI(Tl) crystal calorimeter are combined to form candidates, which must have 
within 2.5cj (about 15 MeV) of m^o. Both daughter photons must be detected in the “barrel” 
region of the detector, have energies of greater than 30 MeV, and deposit most of their energy in 
a compact group of crystals . 

Since D*~^ fragmentation is relatively hard |jl^ and combinatoric background comes mostly 
from low-momentum tracks, we impose a cut of X = p{D*~^)/p(D^g^^) > 0.55. For each event we 
calculate AM, the difference between the reconstructed D*~^ and D'^ masses. We require AM to 
be within 2.5 MeV (3 cj) of the known mass difference. 

Events in which a random slow vr® is combined with a correctly reconstructed D~^ will contribute 
to the peak in M{D~^) |^, but will not peak in the AM distribution. In order to remove this 
background, we perform a sideband subtraction in AM. The resulting invariant-mass distributions 
for all events passing the cuts is shown in Figure 1. 

The reconstruction efficiencies for the signal and normalization modes were estimated using a 
GEANT-based Monte Carlo simulation |]l^ of the CLEO II detector. Furthermore, to study the 
combinatoric background in M{D^) for each mode, we ran a full Monte Carlo simulation that 
included all particle decay processes except for the signal mode and peaking backgrounds from 
other specific decay modes. In all three decay modes, the combinatoric backgrounds are smooth 
and are fit well using a quadratic polynomial. 

Figure 1(a) shows the KsK^ invariant mass spectrum. The peak at about 1.95 GeV is from 
—> Ks't^'^ events in which the vr"*" is misidentified as a kaon. The broad peak in the low mass 
region is from —> Ksp'^, —> Ksp^, and —> K*{892)~tt~^ events. In each of these events, 

a charged pion is identified as a kaon, and the other pion is undetected. The shapes of these 
peaks were obtained from Monte Carlo. The relative normalization of each mode was fixed to 
PDG |10| values, and the overall normalization of the sum was allowed to float in the fit. The 
combinatoric background is parametrized by a quadratic polynomial. The signal is fit with a sum 
of two Gaussians. The ratios of the widths and areas of the two Gaussians are obtained from 
signal Monte Carlo, and the overall width is allowed to float. We find a signal yield of 70.3 ± 12.1 
events at the mass. The reconstruction efficiency, e, is (6.91 ± 0.23)%. As a cross-check, we 
obtain B{D~^ —> KsK~^)/B{D~^ Kstt'^) = 0.28 ± 0.07 from the normalization of the reflection 

background component, which is consistent with the direct measurement. 
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M(Kg7r'^) (GeV) 

FIG. 1. The experimental data and fits. Data are the solid points. The solid lines are fits to 
the data; the dashed lines are the combinatoric background components of the fit. (a) Upper plot: 
M{KsK~^). (b) Lower plot: M{KsTr~^). 
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Figure 1(b) shows the Kst^~^ invariant mass spectrum. The background below 1.75 GeV is 
primarily D~^ which is small and far from the signal, so we exclude this region. The 

region between 1.75 GeV and 1.80 GeV is enhanced by D~^ K^K~^] we obtain the shape of this 

background with Monte Garlo and include it in the ht, allowing the normalization to float. The 
combinatoric background and signal are fitted using the same procedure as above. We observe 
473 ± 26 events with e = (9.32 ± 0.27)%. 

In the normalization mode of D'^ K~7r~^TT~^ we observe 5430 ± 108 events with e = (12.43 ± 

0.19)%. 

The systematic errors are summarized in Table I. To study the particle identification cuts, we 
use a kinematically identified sample of kaons from the decay chain ^ D^Tr~^, K~n^. 

The cut efficiency as a function of kaon momentum for both data and Monte Garlo is measured, then 
integrated over the momentum distribution of Monte Carlo KsK~^ events. This yields 

an overall momentum-weighted efficiency. We find eMc{D~^ KgK^)/edata{D~^ KgK^) = 
1.100 ± 0.030, so a correction factor of 1.10 is applied to the efficiency-corrected yield, N, of 
KsK^. From a similar study of the secondary vertex requirements, we obtain correction 
factors of 1.030 ± 0.014 for N{D^ Kstt^) and 1.036 ± 0.011 for N{D^ KsK^). 

The systematic error is due to differences between the Monte Carlo simulation and 

data in the Dalitz plot distribution of events. The systematic error in the fitting 

procedure was estimated by varying the Monte Carlo background shapes, fitting functions, fit 
regions, and bin sizes. The systematic error for Monte Carlo tracking efficiency is small because 
we measure ratios of branching fractions, and all decay modes have a hnal state of three charged 
tracks. 


Systematic bias 

B(KsK+) 

B(KsTr+) 

B(Ksn+) 
B{K~7T^7r^ ) 

particle ID 

4.0% 

1.0% 

Ks detection efficiency 

1.5% 

2.0% 

Dalitz structure 

0.0% 

3.5% 

7r° cuts 

1.2% 

1.2% 

Sideband determination 

2.7% 

2.7% 

Fitting Procedure 

6.4% 

1.5% 

Signal shape 

2.8% 

2.8% 

Tracking efficiency 

1.0% 

1.0% 

TOTAL 

8.7% 

6.1% 


TABLE I. Systematic error summary. For each ratio measurement, we list each contribution 
to the systematic error in terms of percentage of the measurement. 


The final results are; 


B{D+ KsK+) 
B{D+ Ks7r+) 


0.222 ± 0.041 ± 0.019 


B{D+ KsK+) 
B{D+ K-TT+TT+) 


0.0386 ± 0.0069 ± 0.0037 
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B{D+ Kstt+) 
B{D+ —> A'“7r+7r+) 


0.174 ±0.012 ±0.011 


where the first error is statistical and the second error is systematic. 

To find absolute branching fractions, the last two results are combined with the PDG value 
B[D'^ —> = (9.1 ± 0.6)% |j^. When converting the measurements of Ks branching 

fractions to branching fractions involving we must take into account the possibility of interfer¬ 
ence between doubly Cabbibo suppressed and favored modes [15|. The amplitudes interfere with 
a magnitude of roughly 2 tan^(0c') cos(/3), where 9c is the Cabbibo angle and (3 is an interference 
phase. Because (5 is unknown we assign an 10% uncertainty to B{D~^ —> K^7r~^). There is no such 
interference in D~^ —> K^K~^. We obtain 


B{D+ K^K+) = (0.70 ± 0.12 ± 0.07 ± 0.05)% 


B(T)+ ^ 77%+) = (3.17 ± 0.21 ± 0.19 ± 0.21 ± 0.32)% 

where the third error is due to uncertainty in the normalization branching fraction and the fourth 
error is due to the possibility of DCSD interference. Our result for B{D'^ —> K^K^)/B{D~^ 
X‘^7r+) is 3.6cr higher than iai}?{9c), consistent with the expectation that destructive interference 
suppresses the D+ —> rate [p^ . 

The amplitudes for the three D KK decays may be decomposed into isospin amplitudes; 

= ~^{Ai ± %o) 

= i=(%i-%o) 


^+0 = y/2Al 

where A+- = {K+K-\H\D^), = {K^K^\H\D^), and %+° = {K+K^\H\D+). We have assumed 

that the Hamiltonian reponsible for these decays has isospin structure |/, Is >= |^)±5 >■ This is 
true for all Standard Model decay processes, except for two; the ss-popping W-exchange diagram 
of —> K^K^, and the 11+ annihilation diagram of 71+ ^ .^^77+. Since these diagrams are 

helicity-suppressed and require ss popping, they are expected to be small. 

From these equations one can express the ratio of isospin amplitudes and the isospin phase 
angle difference in terms of measured branching fractions; 

Hi 2 _ r+o 

- 2r+- ± 2roo - r+o 

pH— _ pOO 

cos(5kk) = — , , = 

\/f^%2r+- ± 2roo - r+o 

for D 7777, where 6kk = ai'g(^i/^o)- 
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In the Cabbibo-favored D Kir system, the Hamiltonian has isospin structure |1,+1 >. The 

isospin decomposition and the equations for I As Mil and Sxn = argfAsMi) in the D Ktt 

2 2 2 2 

system are also similar to those oi D ^ KK and may be found elsewhere p7| . 

CLEO has now measured the six branching fractions necessary to calculate the amplitude ratios 
and phase shifts in D ^ KK and D ^ Ktt |§-0. All branching fractions are written in terms of a 
fraction of B = B{D^ —> K~tt~^), in order to avoid additional statistical error from the uncertainty 
in B. We use the CLEO result B{D'^ —> K~tt^tt^) = (2.35 ± 0.16 ± 0.16)B |]^ and the PDG fit 
result B{D^ —> K^TT~^7r~) = (1.41 ± 0.11)B |]T^. The results are listed in Table II. 


Measurement 




Em 

p+0 


D KK 


D ^ Ktt 


(0.116 ±0.010)B 
(0.014 ± 0.004) B 
(0.182 ±0.041)B 


B 

(0.620 ±0.126)B 
(0.819 ±0.136)B 


amplitude ratio 
cos 5 


= o.6i±[]:li 

0.88±[]:JO 


1 

Ai 

T 


= 0.27 ±0.03 


_n 1 2+0.23 


TABLE 11. Isospin analysis inputs and results. r+ = |A+ p, T'^^ = = |A+°p, and 

B = B{D^ K-TT+). 


In conclusion, we find that the isospin phase shift difference in D 
than that of both D — > Ktt and D 
Ai 


KK is significantly smaller 
Eurthermore, the ratio of —> KK 


isospin amplitudes 


a;; 


tttt [5^^ = 0.14T0.16 IH ) 

is 3.5cj from one, which is the value obtained if only Standard Model 
spectator and penguin diagrams are allowed, exchange diagrams are neglected, and the ESI are 
entirely elastic. Therefore, the substantial rate observed for —> K^K^ can be attributed to one 
or more of the following: (1) inelastic rescattering (ESI) from other modes (e.g. —> (pcj) —> K^K^ 

|r^ ) (2) large contributions from W-exchange with ss-popping in decay or (3) annihilation 
diagrams (unlikely in a factorization model [pA[] ). 
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